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Abstract 



There are known models of spherical gravitational collapse in which 
the collapse ends in a naked shell-focusing singularity for some initial 
data. If a massless scalar field is quantized on the classical background 
provided by such a star, it is found that the outgoing quantum flux of 
the scalar field diverges in the approach to the Cauchy horizon. We 
argue that the semiclassical approximation (i.e. quantum field theory 
on a classical curved background) used in these analyses ceases to be 
valid about one Planck time before the epoch of naked singularity 
formation, because by then the curvature in the central region of the 
star reaches Planck scale. It is shown that during the epoch in which 
the semiclassical approximation is valid, the total emitted energy is 
about one Planck unit, and is not divergent. We also argue that back 
reaction in this model does not become important so long as gravity 
can be treated classically. It follows that the further evolution of the 
star will be determined by quantum gravitational effects, and without 
invoking quantum gravity it is not possible to say whether the star 
radiates away on a short time scale or settles down into a black hole 
state. 
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1 Introduction 



If a spherical star collapses to form a naked shell-focusing singularity, the 
nature and magnitude of the energy emission accompanying such a collapse 
is of interest. This emission may broadly be classified into three categories 
- classical |1] (where both matter and gravitational degrees of freedom are 
treated classically), semiclassical (where only the matter degrees of freedom 
are quantized on a classical curved background), and quantum gravitational 
(where both matter and gravity degrees of freedom are quantized). 

A priori, it is not possible to say which of these three kinds of emission will 
dominate. In the present paper we bring together and analyze some of the 
results concerning the semiclassical emission. It is known that the quantum 
Hawking radiation is divergent as the Cauchy horizon is approached. We 
provide a fresh analysis of the divergence of the quantum flux on the Cauchy 
horizon, and derive the emission spectrum for a massless scalar field on the 
spherical dust background. The divergence of the quantum flux has been 
observed in various models [@], ||, ||, || and the naked singularity spectrum 
was studied by Vaz and Witten and by Harada et. al. ||. 

The principal aim of the present paper is to show that during the epoch 
in which the semi-classical analysis is valid, the total emitted energy is about 
one Planck unit and is not divergent. If follows that, during this epoch, the 
back reaction is also unimportant for an astrophysical star because the total 
radiated energy is small compared with the mass of the star. Now the semi- 
classical approximation ceases to be valid about one Planck time before the 
formation of a Cauchy horizon because by then the curvature in the central 
region of the star reaches the Planck scale. Therefore further evolution of 
the star, i.e., during the last Planck time before the formation of the Cauchy 
horizon, will be determined by quantum gravity. An important consequence 
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is that a star collapsing to a naked singularity enters the quantum gravita- 
tional phase with most of its mass intact. Because only about a Planck mass 
has been radiated away until then, semi-classical techniques cannot indicate 
whether the star radiates away a significant portion of its total mass on a 
short time scale or settles down into a black hole state. 



2 Radiated Power 



When a massless scalar field is quantized on the background of a collapsing 
spherical star, the radiated power can be calculated in the geometric optics 
approximation. The expression for the radiated power is J| 
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where P and P refer to a minimally coupled scalar field, 0, and a conformally 
coupled scalar field, 0, respectively. The difference is due to the nature of 
the coupling to the gravitational field and is seen to be only a total diver- 
gence which does not contribute to the total radiated energy under some 
physically reasonable assumptions. In this sense, it may be assumed that 
it is the radiated power for the conformally coupled scalar field that is rel- 
evant for the observed energy flux because this flux is positive definite. P 
and P are the semiclassical contribution to the expectation value of the off- 
diagonal component of the stress-energy tensor, (Tuu)- We note that it is 
entirely expressed in terms of the function G(U), which maps ingoing null 
rays, V =constant, from X~, to outgoing null rays U =constant on X + , i.e., 
V = G(U). This description of the radiated energy in terms of null rays is a 
consequence of the fact that the scalar field modes on the null infinites are 
most naturally expressed in terms of the null coordinates, U and V there M. 



The aforementioned divergence in the radiated power is observed when the 
expressions (|l|) or @ are evaluated for the case of a spherical star collapsing 
to a naked singularity. No such divergence is observed for a star collapsing to 
a black hole. In the following section, we display Q(U) for both cases. In the 
former the function is polynomial, in the latter it is exponential. This key 
difference leads not only to the divergent flux but to a non-thermal spectrum 
as the star collapses to form a naked singularity. 



3 The Dust Solution and the Map Q(U) 



The spherical collapse of dust, as described by the Tolman-Bondi model, 
leads to the formation of a naked singularity for some initial data @, ||, 
[p]|], JTTJ . One can quantize a massless scalar field on the background of such 
a spacetime, and calculate the quantum power radiated in the low angular 
momentum modes, using the geometrical optics approximation described in 
the previous section. This has been done for the self-similar dust model Q], 
for non-self-similar dust collapse [12] and for a C°° model ||, all of which 



are globally naked models. There is some indication from lower dimensional 
models that the locally naked models do not radiate catastrophically fT3|j . 

The marginally bound Tolman-Bondi metric is given by 

ds 2 = dt 2 - R ,2 dr 2 - R 2 dn 2 . (3) 

The evolution of the area radius is given by 



R 3/2 = r 3/2 - ls]F{r% (4) 

where F(r) is the mass function, and the density evolution is determined by 
the equation (we set G — c — 1) 

8^,0 = ^. (5) 
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The cloud is assumed to extend up to some comoving radius r& and is 
matched to a Schwarzschild exterior. One can construct null coordinates 
(u, v) in the interior and the Eddington-Finkelstein null coordinates (U, V) 
in the exterior region. 

The map Q{U) defined in the previous section can be calculated, once 
an interior model for the collapsing spherical dust star is given. Let us first 
consider the self-similar model M. In this case, the mass function is given 
by F(r) = Xr and we have the scaling R — r at the singular epoch, t = 0. 
The density at the center evolves with time as 

becoming infinite at the singular epoch t = 0. 

At the boundary of the star, the Schwarzschild coordinates T, R are re- 
lated to the Tolman-Bondi coordinates by 



fR~ b + y/2M 

3V2M' b -v— -° 1 — v /^r_ V /2M I 
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whence the Eddington-Finkelstein null coordinates become 

Th 4 8 

U b (y) = y 3 - -ar b y - r b y 2 - -a 2 r b In \3y/2a - 1|, 

a 3 9 

V b (y) = - -y 3 - -ar b y + r b y 2 + -a 2 r b In \3y/2a + 1|. (8) 
a 3 9 



where y = y R/r and a = 3\f\/2. 

Tracing a ray V = constant from X~ through the cloud and out again 
toward X + , one finds that, for the self-similar model forming a naked singu- 
larity, the map Q(U) is given by [Q 

Q(U)=A-B(U -U)$ (9) 
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where 7 is a positive constant, less than unity, which is determined by the 
only free parameter, A, in the self-similar model. Above, Uo is the Cauchy 
horizon. On the other hand, if the self-similar star collapses to form a black 
hole, the map is given by fl4 



g(U) = A-Bexp [-U/AM] , (10) 
where M = F(r^)/2 is the total mass of the star. 



4 Spectrum of Created Particles 

Let us evaluate the expectation value of the number operator of the created 
particles, in the Minkowski in-vacuum. This is given by 

roc 

(0|JV( W )|0>= / dJ \(3(u;', u)\ 2 (11) 
J 

with (3 is the Bogoliubov coefficient given by 

P{«>'u) = ±M r dVe-^e-^M. (12) 

Z7T V OJ J-00 

Here, the function U = J~(V) is the inverse of the map V = Q(U). The 
ingoing ray Vq is the last ray which becomes an outgoing ray after passing 
through the center of the star. In case of the naked singularity model, one 
needs to make a subtle assumption about the nature of emission from beyond 
the Cauchy horizon. The present analysis amounts to the assumption that 
there are no outgoing modes from beyond the Cauchy horizon M. In prin- 
ciple, one could consider alternate boundary conditions on and in the future 
of the Cauchy horizon (see for instance Harada et al. f|). 

It follows that, 

(0|iV»|0> = -V r dV [ V0 dV r dJu/e^V'-^e^W-™ (13) 

AtT 2 UJ J -00 J~oo Jo 
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and carrying out the integral over u' we find 

(N( U )) = f° dV f V ° dV> - * e ^0-W]. (14) 

4tc z u J~oo J-oo (V ' — V + ie) z 

Now inverting the expression (|T0| ) for Q(U) to obtain the map FiV) for the 
black hole, and making the substitutions x = V$ — V, y = Vo — V, and then 
p = x — y, q = y/x we get 

1 r°° dp r°° q-toMu 
v ' 4ttW-oc p Jo (1-q + ie) 2 



The integral over q can be performed [HJ. Finally taking the e — > limit, 
we obtain 



iVH = — ^r-r/ -A (16) 



2M 1 r 00 rfp 
7T e 87rM£J - 1 7-oo p 

Thus we recover the expected thermal result for the frequency dependence 
in the case of the black hole. The divergence resulting from the integral over 
p represents a steady flux over an infinite time and can be eliminated by 
considering the emission per unit time. 

The integral for the naked singularity case can be performed by inverting 
the map given in @ and by making the substitutions x = V — V, y = V — V 
followed by p = x 7 — y 7 , q = y/x. This gives 

e iu>Bp j 



N(u) - 1 



roo poo rU ri 

/ dp dq + dp dq 
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e- 1 
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or, after integrating, 

The frequency dependence of the spectrum is 1/u. In evaluating the integral 
( |14D for the naked case we have made the change e — > ex. As is seen, the role 
of the regulator e in the naked case is quite different from that in the black 
hole case. The limit e — > yields a divergence - the origin of this divergence 
appears to be different from the black hole case, since now the flux is not for 
an infinite time. 
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The fact that the divergent behaviour of the spectrum does not depend 
on the parameter 7 characterizing the naked singularity suggests that this 
spectrum may be the result of imposing self-similarity on the solution. It 
could also be related to the cutoff at the Cauchy horizon (i.e. the assump- 
tion that there are no outgoing modes beyond the Cauchy horizon). It also 
suggests that an inference of the spectrum will strongly depend on the nature 
of boundary conditions imposed on the Cauchy horizon. 

We would like to emphasize that the boundary conditions assumed here 
for the spectrum calculation are not at all relevant for the calculation of the 
radiated energy described in the next section. The result on the radiated 
energy is entirely based on calculations in the spacetime region prior to the 
Cauchy horizon. Hence there is no need to make any assumptions regarding 
boundary conditions on the Cauchy horizon, in so far as the calculation of 
emitted energy is concerned. 



5 Radiated Energy 



Let us now calculate the radiated energy. Using (|l|) and the forms (|10|) and 
(Bp for G(U) we get the radiated power in the black hole case as 

h 



p = p 



and in the naked case as 
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The net, integrated, flux emitted in the naked singularity model can be 
written as 
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As seen above, in the limit that U approaches the Cauchy horizon Uq, 
the emitted flux diverges. However, this divergence needs to be interpreted 
carefully. We are working with the semiclassical approximation, in which the 
gravitational field has not been quantized. The semiclassical approximation 
ceases to be valid if the curvature at some point inside the star approaches 
Planck scale. From (|^) we see that this happens at the center of the star 
at t ~ —t pi anc k where £p; anc fc is the Planck time. Hence, we may say that 
the results of the semiclassical approximation cannot be trusted beyond this 
epoch. 

From equations (^) and (§) one can show that for a ray starting from the 
origin at t — —tpi anc k, the time difference Uq — U is greater than of the order 
of Planck time. This is done as follows. Let ray A start from the origin r = 
at t = —tpi and ray B start from the origin at the singular epoch t = 0. The 
outward propagation of a radial null ray is described by the equation 



Since ray B arrives at a comoving coordinate r later than ray A, it follows 
from this equation that the value of dt/dr at a given r will be higher for B 
compared to A. Hence the time interval At = tp — t& between the arrival 
times tp and of these two rays at the boundary of the star will exceed 



Using equation (^[) and the definition of y we can write the difference 
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AU = U Q - U = Up - U A as 



AU = tp — tA + other positive terms 



(25) 
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This shows that the asymptotic time difference Uq — U will be greater than of 
the order of Planck time. Hence, one can see from that in the range that 
the semiclassical approximation can be trusted, the emitted semiclassical flux 
is less than of the order of Planck energy, which is in general much smaller 
than the mass of the original star. 

This suggests that very early during the collapse, the star enters the 
quantum gravitational phase, with most of its mass intact. The further 
evolution, including the nature, duration and amount of emission, will be 
determined by quantum gravity. In this sense, the naked singularity system 
behaves very differently from the black hole system. For an evaporating black 
hole, quantum gravitational effects become important only during the very 
final stages, when the mass of the black hole approaches Planck mass. For an 
astrophysical black hole, the time taken to reach this quantum gravitational 
phase is enormous (1O 71 (M/M ) 3 s). Also, it is an open issue as to whether 
the black hole evaporates entirely, or leaves behind a Planck mass remnant. 

On the other hand, for the star forming a naked singularity, quantum 
gravity becomes important much before the completion of the collapse. And 
it is at present an open issue as to whether the naked star settles into a black 
hole or explodes on a very short time scale. This issue has to be decided by 
quantum gravity. Perhaps one can say that for the first time physicists have 
encountered a dynamical system whose evolution, within the lifetime of the 
Universe, can only be understood with the help of quantum gravity. Indeed, 
the Universe may well contain such stars, whose final evolution cannot be 
understood without quantum gravity. A study of such systems appears to be 
a useful probe of theories of quantum gravity like canonical general relativity 
and string theory. 
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6 Back Reaction 



An important aspect of the above analysis is that we have ignored the back- 
reaction, and one could well ask if our conclusions could be affected by its 
inclusion. We now give arguments as to why the back reaction cannot be 
important during the semiclassical evolution, in the present case. One can 
propose two different criteria for deciding as to when the back reaction be- 
comes important. One is that the total flux received at infinity becomes 
comparable to the mass of the collapsing star. As we have seen above, if the 
mass of the collapsing star is much greater than the Planck mass (as is of 
course usually the case) then the back reaction does not become important 
during the semiclassical phase, because the net emission cannot be more than 
about one Planck mass. 

The second criterion for the back-reaction to become important is that the 
energy density of the quantum field becomes comparable to the background 
energy density, inside the star. In a four- dimensional model, it is difficult to 
establish whether this happens. However, one can study the evolution of the 
quantum field inside the star using the 2-d model obtained by suppressing 
the angular part of the 4-d Tolman-Bondi model |Tj|. Using this model, 



Iguchi and Harada [O have shown that during the domain of validity of the 



semiclassical approximation the expectation value of the energy momentum 
tensor of the scalar field remains smaller than the energy density of the 
dust matter in the collapsing background. Thus the back reaction does not 
become significant during the semiclassical evolution, inside the star, for this 
2-d model. Since a solution of the four dimensional model obtained using 
the geometric optics approximation is equivalent to the exact solution of the 
two-dimensional version, it is plausible that this result holds for the four 
dimensional stellar dust model as well. In this case we can conclude that 
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quantum gravitational effects are more important than the semiclassical back 
reaction in deciding the evolution of the star. 

Again we see that the semiclassical evolution is very different from the 
black hole case, since the back reaction plays a crucial role during the semi- 
classical history of the black hole. 

7 The Smooth Case 

Similar conclusions about the nature of the evolution hold also for a C°° 
model studied by Harada et al. 0]. They also find that the amount of energy 
emitted during the semiclassical phase is of the order of Planck energy. Hence 
our results are not specific to the self-similar model. 

The situation is similar but more complicated than the self-similar case. 
The C°° case is given by 

F(r) = F 3 r 3 + F 5 r 5 + ■■■. (26) 

Then Harada et al. found that the function Q behaves as 

Q = V - A(U -U)- fvl /2 (U - Uf /2 + higher order terms, (27) 

where u s is determined by the initial density profile and the above formula 
is valid only for z/ s (£/q — U) < 1. Using the above formula, they also found 
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It is noted that v s may exceed the Planck frequency vpianck because v s is 
determined by some combination of two macroscopic scales, the free-fall time 
and the scale of inhomogeneity. Nevertheless, we can find that the total 
radiated energy for this late time radiation cannot exceed the Planck energy 
fovpianck ~ 10 16 erg before the Planck time t P i anck = lCT 43 s. The proof is the 
following. If u s > upi anc k, then it is found that the late time radiation is only 
for U Q - U < v~ l < tpi anck . Therefore, we only consider v 8 < vp Xanck . For 
this case, it is found 

Eu -u~t Planck ~ h(y s vpi anck ) l l 2 < hupi anck , (33) 
Eu ~u~tp lanck ~ hv Planck ( — In Planck \ <hvpi anck) (34) 

V ^ Planck v s / 

where the second inequality comes from the fact that the function — xlnx 
has a maximum value e -1 in the domain < x < 1. 



8 Discussion 



One could conjecture that this picture (i.e. the semiclassical flux is of the 
order of Planck energy) will in general hold for a naked singularity forming 
in spherical collapse, even if the matter model is other than dust. This is 
because the curvature induced quantum particle production is a kinemati- 
cal effect and the only role dust plays is that of providing the background 
gravitational field. Any other matter model giving rise to a spherical naked 
singularity can be expected to have a similar semiclassical evolution. The 
fact that the semiclassical emission is so small appears to suggest that the 
emission comes from a very small central spacetime region in the neighbor- 
hood of the naked singularity. If so, then this is a property of the spherical 
geometry, which should continue to hold when pressure is included. In the 
case of a black hole, the semiclassical emission is a property of the event hori- 
zon - as a result the evaporation rate is slower and the semiclassical phase 
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lasts for a much longer duration. 



In contrast with the spherical case, if the collapse is highly non-spherical, 
say resulting in a spindle naked singularity, it is possible that the semiclas- 
sical emission could be a significant fraction of the initial mass. This is an 
important open problem which deserves to be investigated. 

Understanding the quantum gravitational evolution of a spherical star 
which terminates in a naked singularity is thus an outstanding challenge for 
candidate theories of quantum gravity. One would like to know whether the 
star settles into a black hole state or instead, undergoes a catastrophic explo- 
sion on a short time scale, resulting in a cosmic burst. A first step towards 



modelling this phase has been taken in [18 where a formalism is developed for 



the midisuperspace quantisation of the Wheeler-deWitt equation describing 
spherical dust collapse. 
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